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Abstract. More than 30 years after the discovery of the J/ip, its production mechanism 
is still poorly understood. With the LHC data it will be possible to study quarkonia up to 
very high transverse momenta and with high statistics. In this note we discuss experimental 
aspects of J/4' and T(1S) production at the LHC. In particular, we investigate the sensitivity 
of a general purpose LHC detector to observables, which are complementary to the cross 
section and polarisation measurement. These observables would be sensitive to the radiation 
produced in association with the quarkonium. 
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1. INTRODUCTION 



^ ■ Although J/tp^s and T's have been studied extensively since their discoveries, 

cn ! their prompt production mechanism is still not well understood. For reviews 

^ I on quarkonium production the reader is referred to Refs [l|, |2|. Below we only 

fr^ ■ summarise some relevant aspects. 

r^ ! Initially, prompt production of J/ip^s and T's was assumed to take place via 

O I the leading order colour singlet model (CSM). However, in 1997, CDF studies 

5l? ■ showed large excesses in the amount of prompt quarkonia produced with respect 

to theory predictions [3|, |j] , and new production mechanisms were invented. Among 
the most popular approaches is Non-Relativistic QCD, where production at parton 
/\ • level can also take place via a colour octet quark pair. While this colour octet mech- 

j^ ■ anism (COM) is quite successful in explaining the transverse momentum spectra 

of quarkonia, the polarisation predictions of the COM mechanism are in sheer dis- 
agreement with experimental data: production via a colour-octet state predicts a 
significant transverse polarisation of J /ip^ s and T's, while recent Tevatron data 
contradict this prediction |5|, l6|. The CSM in its own has been revived recently, 
when it has been shown that NLO and NNLO QCD contributions [7|,|8|,|9| as well as 
s-channel cut contributions |lOl . Ill| arising from relativistic corrections could lead 



to an enhancement of the production cross section of the CSM, while predicting a 
polarisation closer to the data. 

Given the puzzling situation, it is the utmost importance to study prompt 
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quarkonia at the LHC, where they can be studied with vast statistics up to high 
transverse momentum. Standard measurements are generally related to kinematical 
distributions of the quarkonium decay products, such as differential cross section 
and polarization measurements, and focus on decays into muons. Although these 
provide useful information, they may be sensitive to other factors. Hence it is 
important to investigate additional observables. In particular, we should take 
into account not only the kinematics of the quarkonium itself, but also that of 
particles produced in association. In fact, in some cases prompt quarkonia can be 
searched for in association with a heavy quark pair [91]. Also, for an unambiguous 
understanding of quarkonium production, direct production should be studied, for 
example Xc production could be studied by radiative decays into J/ip^s. 

In this note we discuss general aspects of prompt J/ip and T(1S) production at 
the LHC. Also, we show examples of observables that are sensitive to the hadronic 
activity directly around the quarkonium, thereby allowing to extract information 
about the radiation emitted off the coloured heavy quark pair during production. 
The results apply only to a general purpose LHC detector, and for this work we 
took the CMS parameters 12] as reference. We do not distinguish here between 



direct and indirect (via radiative Xc decay) prompt J/'j/'-production. 

The outline of this note is as follows. In Section 2 we discuss event generation 
of J/ip and T(1S) in PYTHIA. In Section 3 we discuss the production of J/ip^s 
at the LHC, and study the differences in radiation between different production 
toy models. We will also study the J/?/^-reconstruction efficiency. In Section 4 we 
concentrate on T's. Section 5 contains conclusions. 



2. QUARKONIUM EVENT GENERATION 



We have produced J/i/j-events using PYTHIA 6.409 [13|]. The implementation \i£ 
is based on the NQRCD approach and quarkonium production via the singlet (LO) 
and octet mechanism. The octet cc-states have a mass of 3.1 GeV. The values of the 
NRQCD matrix elements are taken from Ref. 15|. We have modified the default 



branching ratios for the radiative decays of the x-states into J/^'s to correspond 
to the recent PDG values [16i]. The polarisation has been set to zero. Production 
via X states is included. In this study we forced the J/ip and T to always decay into 
H~^IJ.~, the branching ratios are 5.93% and 2.48% respectively [l6|]. The following 
features are included. 



2.1. Cross section regular isat ion 

We have made use of an event reweighting function, allowing the quarkonia 
cross section to be dampened at small Pt values, in analogy to the underlying 
event formalism in PYTHIA [l7|, [l8| . The events obtain a weight Wi according to 
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where pT is the transverse momentum evaluated in the rest frame of the scattering, 
Q^ the momentum transfer scale, and pjb the regularisation scale, used also in the 
underlying event cross section regularisation. As can be seen from Eq. [H at small 

Prp a suppression is obtained while at large Prp the cross section is unmodified. 
The value of Ptq is be extrapolated from the Ptq value determined at 1960 GeV 
to the LHC energy: 

/ E \PARP(90) / l4TeV \0-i6 

(2) 

In Fig. [1] (left) we demonstrate the regularisation procedure, as well as the difference 
between the old PYTHIA (before version 6.324) and the new production model. 



2.2. Shower activity associated with quarkonium 

production 

Below we discuss J/'i/'-production, however it applies to T's too. If the J/ip is 
produced via the singlet model, the cc does not radiate and the charmonium 
is thus produced in isolation (when produced at LO), up to normal underlying 
event activity. When it is produced via the octet mechanism, the situation is not 
apriori clear, but it is reasonable to assume that a shower can evolve and that 
the charmonium is embedded in a certain amount of shower activity. This results 
from 1) repeated g —* gg branchings before the branching g -^ cc, and 2) from 
the radiation off the colour octet cc-state once formed. There are several switches 
to regulate the hadronic activity. Using these, we generate J/ip^s in four different 
production mechanisms, each with varying amount of radiation. 

• Production via singlet only. The J/ip^s are produced in isolation. In the hard 
process a number of hard gluons is emitted for colour neutralization, but these 
are recoiling against the J/ip and thus expected to be back-to-back to the J/V'. 

• Octet production only, where the octet cc undergoes only a small amount of 
parton shower evolution only. Here MSTP(148)=0, implying that the splitting 
kernel used for the branching cc*^^^ -^ cc^^^g is the q ^ qg kernel, modified for 
the cc-mass. Only a small amount of activity if expected to be emitted during 
the parton cascade, yet the activity is expected to be much larger than in the 
colour singlet case. This case is referred to as octet low radiation scenario. 

• Octet production only, but MSTP(148)=1 and MSTP(149)=0. Now the split- 
ting kernel for cc*^^-* -^ cc^^^g is the g ^ gg one, again modified for the cc-mass. 
In each branching the cc-state is likely to get the largest energy: z > 0.5. The 
amount of radiation is larger than in the previous case, however because of the 
suppression of z < 0.5 only a small amount of hard gluons is expected around 
the j/ijj. This case is referred to as octet medium radiation. 

• Octet production only, but now MSTP(148)=1 and MSTP(149)=1. Again 
the g ^ gg splitting kernel is used, but now z > 0.5 and z < 0.5 have equal 
probability (in reality z > 0.5 still dominates due to mass effects). Here a 



significant amount of liard gluons can be emitted before the hadronisation 
process, in fact showering activity might be overestimated. We will refer to 
this case as octet high radiation scenario. 

It must be emphasised that these models are 'straw- man models', and that we do 
not expect any of them to reproduce the data by itself. Finally, note that CDF 
data can be reproduced by adding singlet and octet contributions, but that any of 
the 3 octet models together with the singlet model fit the data. 



2.3. Event generation parameters 

For the production of prompt J/^— > fi^fi^ -events we have generated events in 5 
bins oi pT, the transverse momentum exchanged in the hard process (0-10, 10-20, 
20-30, 30-50, 50-inf GeV). This was done in view of the large production cross 
section at low Prp , in order to have statistics also at large Prp' . In total 2M 
events in four different production scenarios (see Sec. 12. 2p have been generated. 
The luminosity of the lowest px-bin for the singlet model corresponds to only 0.2 
pb~^, while the higher pxbins correspond to higher luminosities. We applied a filter 
at generator level requiring two muons with \ri\ < 2.5 and pt > 2 GeV/c. 

As background for the J/'?/'-production studies we have considered non-prompt 
J/i/j-events. An inclusive non-prompt J/tp sample has been generated (PYTHIA 
setting MSEL=1, cross section 55 mb), and passed through a dimuon generator 
level filter (|?7| < 2.5 and pt > 2.5 GeV/c). The luminosity of the sample we use is 
4.45 pb"^, corresponding to roughly 400000 filtered non-prompt J/psi-events. 

Prompt T's have been generated in the four different production scenarios of 
Sec. 12.21 corresponding to 100 pb~^, without pT-bins. Here a dimuon filter was used 
requiring the muons to have |?7| < 2.5 and pt > 3 GeV/c. Since the T-production 
cross section is smaller than that of J/ip^s, we will need more LHC data, possibly 
including pile-up. We included here 5 pile-up events per T-event. 

As background for the T-production studies we have used an inclusive pp — * /iX 
sample generated again with MSEL=1 in PYTHIA, corresponding to 0.14 pb^^. 

All events are processed through a toy-detector based on the expected perfor- 
mances of the CMS detector |12l]. We have not included trigger requirements, but 
being interested in high values of the transverse momentum of the quarkonium (see 
Sec 4), we expect no problems in triggering. Typically, a double muon with some 



threshold would be used, for example 3 GeV for each muon [12 



3. J/ij PRODUCTION STUDIES 
3.1. J/i/j production at the LHC 

In Fig. [1] (right) the differential cross section for a CMS-like detector as function 
of Prp'^ is shown for the four production cases. To obtain the total cross section the 
singlet curve should be added with one of the three octet curves. The singlet curve 
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FIGURE 1. Left: J/ip's at the Tevatron. The old model (singlet only, normahsed to fit CDF 
data, green line), the new model without reweighting (blue broken line), the new model with 
reweighting (red line), and the CDF data (black data points, from Ref. 1^). Right: differential 
cross section at CMS-like detector for the four production cases separately. 



decreases rapidly with Prp' , as expected. For the octet cases, the decrease is faster 
for models with more radiation, since here the cc loses more energy in radiating. 
Although the curves are well separated at high transverse momentum, these curves 
are influenced by different factors, such as initial and final state radiation, the mass 
of the cc-octet, and parameters entering in the cross section dampening [2Q[. In 
conclusion, it is important to find complementary observables. 

To have a rough idea of the statistics expected at CMS, we display in Tab. [1] the 
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numbers of J/ip^s produced within acceptance in ranges of Prp' in 100 pb 

Concerning the production of non-prompt J/ip^s, the rate at LHC is significant, 
as can be seen in Fig. [2l While at low Pr,/ prompt production dominates, at higher 
Prp' the non-prompt and prompt contributions are comparable. 



TABLE 1. Rough estimate for number of prompt 
J/V''s produced and reconstructed in 100 pb~^ (singlet 
and octet medium radiation summed). Reconstruction of 
J/ip's is discussed in Sec. 



P^/'^-range (GeV) 


Nr. produced 


Nr. reconstructed 


0-10 


3.0x10^ 


3.7x106 


10-20 


3.5 X 10^ 


1.6 xlO^ 


20-30 


2.0 xlO^ 


1.4 xlO^ 


30-40 


3.2 X 10^ 


2.2 xlO^ 


40-50 


7.8 X 10^ 


5.6x10^ 


>50 


3.5 X 10^ 


2.4x10^ 



Prompt J/\|/: singlet + octet medium radiation 
Non-prompt J/y 
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FIGURE 2. Differential cross section for non-prompt (black) and prompt (red) J/i/i's. 



3.2. J /ifj reconstruction 



Reconstruction of J/V'-events with the CMS detector is discussed in Ref. [21 



We display in Fig. [3] the reconstruction efficiency with respect to the filtered events 
of our toy detector as function of Prp . As long as no isolation requirements are 
applied, the reconstruction efficiency is independent of the production model. It 
can be seen that J/V^'s with Pt below 3 GeV are hardly reconstructable in a 
detector like ATLAS or CMS, because these J/ip^s decay into two approximately 
back-to-back muons with too low momentum. It must be noted that both ALICE 
and LHCb are able to reconstruct J/ip^s with much lower transverse momentum. 
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FIGURE 3. Left: expected J/tp reconstruction efficiency in general purpose LHC detector as 
function of P^ . Right: the same but for T's. (see Sec. l4.2p 



3.3. Sensitivity to alternative J/^/'-production observables 



Since the shower-activity is different for the different models, we investigated 
several variables sensitive to the kinematics of the particles produced in association 
with the J/ip- Because at lower Prp , shower activity is small in general, differences 

manifest themselves only at higher values of P-p , about 20 GeV/c. It must be 
noted that the energy of the surrounding particles associated with J/-?/'-production 



is small (order GeV), and as such it is not apriori clear whether there is any 
sensitivity at all. 

Examples of observables were shown in Ref. 22]. In Fig. H] (left) we dis- 
play for the four prompt J/psi production models the fragmentation variable 



pJ/V>/^pJ/^ 



+ PS~°-'^) as function of P^ , after reconstruction criteria. Here 



pR=0.7 



is the sum of the transverse momentum of particles in a cone oi R = 0.7. 
We remark that z is influenced quite strongly by the underlying event activity. In 
fact, without underlying event activity z f^l for the singlet case. 

In Fig. m (right) we display the transverse momentum density dPr/dfi^ for J/ip^s 
between 20 and 40 GeV after reconstruction, in a cone around the J/tp of certain 
size R = ^J{^r]Y 



(A0)^, where 



J p around 
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paround 



(R + dR/2) 
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rrp 



(R-dR/2) 



d(^R 



7r[(R + dR/2)2-(R-dR/2)2 



(3) 



Here P^^°^^^(R) is the sum of the transverse momentum of all charged particles 
(with Pt > 0.9 GeV) inside the cone of size R. At this stage, we can conclude that 
the particles produced in association with prompt J/^-production, even though 
low energetic, are well detectable with a general purpose detector at LHC. 




J/psi activity: 20 GeV < P"'* < 40 Ge^ 
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FIGURE 4. Left: variable zj/.^ at reconstruction level. Right: variable — ^^ — —^ at recon- 
struction level. 



3.4. Background 



The main backgrounds for prompt J/'i/'-production are non-prompt J/ip^s and 
fake J/ip^s from QCD background. The latter can most probably be easily sub- 
tracted with sideband subtraction, but for non-prompt J/^'s the situation is tricky. 
By studying the impact parameter of the reconstructed J/ip^s, and knowing the 
behaviour of prompt and non-prompt J/ip^s, it is possible to evaluate the amount 
of non-prompt background in data. However, to correctly subtract the non-prompt 
contribution from data, one must know very precisely the hadronic activity around 
non-prompt J/z/^'s. It turns out that the procedure is difficult, in particular be- 
cause the hadronic activity around non-prompt J/ip^s is very large in comparison 
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FIGURE 5. Left: T production with singlet and octet in PYTHIA (red), together with CDF 
data (back dots). The latter are taken from Ref. ^^. Right: Differential cross section at CMS-like 
detector for the four production cases separately. 

TABLE 2. Estimate for number of prompt T's pro- 
duced and reconstructed (singlet and medium radiation 
model added) in 100 pb~^. Reconstruction of T's is dis- 
cussed in Sec. [ 



P^^'^-range (GeV) 


Nr. produced 


Nr. reconstructed 


0-10 


1.3x10^ 


5.7x10^ 


10-20 


2.1 X 10^ 


1.4x10^ 


20-30 


3.0 X 10^ 


2.3 xlO* 


30-40 


6.1 X 10^ 


4.8 xlO^ 


40-50 


1.7x103 


1.4x102 


>-50 


9.1 X 10^ 


7.0x102 



with that of prompt J/ip^s. When attempting to subtract the background, the 
differences between the prompt J/7/^-production models were washed out. Several 
improvements are possible, such as more Monte Carlo statistics, or by combin- 
ing different variables, to gain sensitivity. However, given the fact that for T's no 
non-prompt background is present, we first investigate hadronic activity there. 

4. UPSILON PRODUCTION STUDIES 
4.1. T production at the LHC 



We generated T(lS)-events as described in Sec. I2.3[ In Fig. [5] (right) the differ- 
ential cross section for prompt T(1S) at LHC as function of P^ is shown for the 
four production cases. In comparison with the J/i/j-ctoss section in Fig. [H we see 
that the curves for T's fall less rapidly than for J/ip^s. Furthermore, although the 
difference between the singlet and octet curves is still large, the differences between 
the three octet models is small, and the curves are almost overlapping. The reason 
is that a coloured bb-state emits less gluon radiation than a cc because it is heavier. 

A rough estimate of the statistics expected at a general purpose detector like 
CMS or ATLAS is given in Tab. El 



4.2. Upsilon reconstruction 

Just like J/ip^s a general purpose detector could detect T's by detecting the 
muons. The reconstruction efficiency with respect to the ffitered events as a function 
of P^ if displayed in Fig. [3] (right). As opposed to J/ip's, T's can be reconstructed 
down to zero transverse momentum. Indeed, even when P^— > 0, the momentum 
of the then back-to-back muons is sufficiency large to be detected, because of the 
high T-mass (9.46 GeV). Furthermore, the value of the reconstruction efficiency 
where it plateaus is seen to be larger for T's than for J/ip^s. This is due to the 
fact that the muons from T — > fi'^fi~ have a larger separation angle than those of 
J/\l' — > //"""/i", so are better detectable. 



4.3. Sensitivity to alternative observables 

In Figl6] (left) we display the variable zy at reconstruction level. In Figj6] (right) 
the transverse momentum density dPr/dfi/j is shown for T's between 20 and 40 
GeV for the four T production models, in a cone R around the T, at reconstruction 
level. Although the amount of hadronic activity is smaller than for J/^'s, the 
differences are still well visible. To obtain these results, ffist the tracks coming 
from pile-up have been subtracted. 
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FIGURE 6. Left: variable z~r at reconstruction level for T's. Right: variable 
reconstruction level for T's. 
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4.4. Background 

The main background is from QCD events, where we have b-quarks decaying 
into muons. The spectrum is well-known. Due to the weak b-decay, the hadronic 
activity around these reconstructed muon-pairs is high. Unfortunately, studying 
this background with the current Monte Carlo samples is difficult. The efficiency 
for app—y fiX background event to end up both close to the T-mass and at high P^, 
is very tiny. To study the background with sufficient statistics, a more dedicated 
Monte Carlo generation should be done. In real data, however, we expect studying 



background not to be a limiting factor to this analysis, because it can be studied 
accurately in the sidebands left and right of the T-mass window. Having no other 
background components, most probably it can simply be subtracted to obtain the 
hadronic activity of pure T's. 

5. CONCLUSIONS 

Up to today the prompt quarkonium production mechanism is not understood. 
LHC will open up new possibilities for studying quarkonium production. Apart 
from cross section and polarisation measurements, it is important to study different 
observables. An example is shown here, and suggests studying observables sensitive 
to hadronic activity produced in association with prompt quarkonia. 
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